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ABSTRACT 

A design procedure for subscale airfoils with full- 
scale leading edges that exhibit full-scale water 
droplet impingement characteristics in an incom- 
pressible, inviscid flow is presented. The design pro- 
cedure uses validated airfoil design, flow analysis and 
water droplet impingement simulation codes to ac- 
complish the task. To identify and isolate important 
design variables in the design, numerous trade stud- 
ies were performed. The paper presents the results 
of the trade studies and briefly discusses the role 
of important design variables in the subscale airfoil 
design. The effect of these design variables on circu- 
lation, velocity distribution and impingement char- 
acteristics is discussed along with the accompanying 
implications and compromises in the design. A strat- 
egy to incorporate viscous effects into the design is 
nlfln presented. The paper also presents the design 
of a h alf -scale airfoil model with a 5% upper and 
20% lower full-scale surface of the Learjet 305 airfoil 
leading-edge and compares its aerodynamic as well 
as the droplet impingement characteristics with that 
of the Learjet 305 airfoil. 

NOMENCLATURE 

c = airfoil chord length 

C 4 = airfoil drag coefficient 

C\ — airfoil lift coefficient 

Cfno = airfoil pitching moment coefficient 

F r = Froude number, Ufy/cg 

K = droplet inertia parameter, p w 6 2 U/\8cp 

Ks = trailing-edge thickness parameter 

M = freestream Mach number 

Re = freestream Reynolds number, pUc/p 

Ry = droplet freestream Reynolds number, 

pSU/p 

S = airfoil surface arc length measured from 

the leading-edge 
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T = freestream static temperature 

u, v = local flowfield horizontal and vertical 

velocity components 
U = freestream velocity 

V = surface velocity 

VM D = volume median droplet diameter 


Xy y = airfoil coordinates 
x 0% y 0 = initial horizontal and vertical 
displacement of the droplet 
x m , x m = upper and lower surface match locations 
Xri x r = upper and lower surface pressure 
recovery locations 

Vx — design velocity level for segment 1 

a = angle of attack relative to the chord line 

a t = effective angle of attack relative to 
the nose section chord line, a — 7 
a * a* = upper and lower surface multipoint 
design angle of attack distribution 
0 = local impingement efficiency 

<f> u = leading-edge arc limit 

y zz nose droop angle 

T = circulation strength normalized by Uc 

Y = circulation strength, m 2 /s 

S = droplet diameter 

rj = normalized subscale airfoil chord 

length, /c/. 


= air viscosity 
= air density 
= water density 
ss finite trailing-edge angle 


Subscripts: 

fs — full-scale airfoil 

i = inviscid 

/ = lower surface 

ss = subscale airfoil 

u = upper surface 

11 = viscous 


INTRODUCTION 

Recent aircraft accidents have raised important 
flight safety issues related to the effect of ice accre- 
tion on airfoil and wing performance. In order to 
improve flight safety, a better understanding of the 
effect of ice accretion on the aerodynamic perfor- 
mance of modern airfoils is required. One important 
step in this process is to evaluate the aerodynamic 
performance of the airfoil sections (or the wing as a 
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whole) at the icing conditions within the certification 
icing envelop that result in the largest performance 
penalties. 

Since ice accretion scaling is still not well under- 
stood, testing at full-scale or near full-scale condi- 
tions is highly desirable. The available ice accretion 
tunnels, however, are too small to test full-scale air- 
foils or wings of most aircraft of interest. Numerous 
investigators have performed experimental or ana- 
lytical studies 1 ” 4 in an effort to evaluate full-scale ic- 
ing protection systems for wing sections using trun- 
cated airfoil models. These truncated airfoil models 
utilize a full-scale leading edge section followed by a 
faired or flapped aft section that, in effect, reduces 
the overall length or chord of the model. To our 
knowledge, however, no systematic study has been 
performed to provide insight into the design of the 
aft section. 

With these issues in mind, a subscale model de- 
sign procedure was formulated with the objective of 
providing design guidance for subscale models that 
simulate full-scale water droplet impingement char- 
acteristics. It is assumed that ice accretion will be 
the same if droplet impingement, surface geometry 
and surface flowfield are the same provided the same 
cloud properties, model surface quality, model sur- 
face thermodynamic characteristics exist. Using the 
fact that ice usually accretes only on the airfoil lead- 
ing edge, where the supercooled water droplets im- 
pinge and form ice, the subscale airfoil model is de- 
signed with a leading-edge geometry (first 10-20% of 
chord) identical to that of the full-scale leading edge. 
The design of the aft section is such that it provides 
full-scale flowfield and droplet impingement on the 
leading edge. Using this formulation, the effect of 
various design variables on the inviscid flowfield and 
droplet impingement characteristics of the subscale 
airfoil was examined to obtain useful guidelines for 
the design. The final design is based on viscous con- 
siderations as well. 

The model design procedure for full-scale flowfield 
and droplet impingement simulation uses validated 
computational airfoil aerodynamics and droplet im- 
pingement codes, 5 ” 15 specifically, an inverse de- 
sign method, 10 the Eppler code, 8,9 XFOIL 11 and 
AIRDROP. 12 

DESIGN APPROACH 

A conceptual illustration of the subscale airfoil de- 
sign procedure is shown in Fig. 1. First, a droplet 
impingement code can be used to predict the lim- 
its of the droplet impingement, which defines the 
initial ice accretion envelop. (The droplet impinge- 
ment code, AIRDROP, 12 is discussed later.) Once 
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Fig. 1 A conceptual illustration of the subscale 
airfoil design procedure. 

the limits of impingement are known over the leading 
edge of the full-scale airfoil, that part of the full-scale 
airfoil geometry is fixed for the subsequent subscale 
airfoil shapes. For the sake of discussion, this fixed 
leading-edge section, which is common to both the 
full scale airfoil and the subscale airfoil, is referred 
to as the nose section while the remaining section of 
the subscale airfoil profile is referred to as the aft sec- 
tion. The aft section of the subscale airfoil is then 
designed to provide full-scale flowfield and droplet 
impingement on the nose section of the subscale air- 
foil. 

An initial geometry for the aft section is obtained 
through the use of a multipoint inverse airfoil design 
code 10 (PROFOIL). The design of this intermediate 
airfoil, from which the aft section of the subscale 
airfoil is derived, is governed by several constraints, 
namely, the scale of the subscale airfoil, the upper 
and lower surface thickness and slope at the junction 
between the nose and aft sections (x m , x m ), and a 
desired form for the pressure recovery characteris- 
tics. Apart from these constraints, additional conti- 
nuity and closure constraints that form an integral 
part of the inverse design methodology 10 are also 
satisfied in order to achieve a physically possible de- 
sign. A multi-dimensional Newton iteration scheme 
is employed to satisfy these constraints. The depen- 
dent and independent Newton variables 10 used in 
the design are listed Table 1. Once the constraints 
are satisfied, the aft section is combined with the 
nose section to form a subscale airfoil. 

The potential flow over both the subscale and 
the full-scale airfoils is then analyzed using the 
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Table 1 Newton variables used in the design. 


Dependent Variables 

Independent Variables 

A'j = 0.3 

<t>le 

c mo 

Vl 

y(*m) 

Q* 

y( x m) 

cT 


Eppler code, which has the capability to analyze 
the potential flow over the airfoils using a method 
that employs panels with distributed surface singu- 
larities. The singularities used are vorticities dis- 
tributed parabolically along each panel. Results pre- 
dicted by the Eppler code have been shown to com- 
pare well with experiments. 16,17 

In order to have a physically similar flow in the 
vicinity of the nose section of both the subscale and 
the full-scale airfoils, the analysis is performed at 
the same angle of attack relative to the nose sec- 
tion chord of both the airfoils. The local inviscid 
velocity distributions over the nose section and the 
stagnation point locations on both the subscale and 
full-scale airfoils are then compared. If the desired 
velocity distribution over the nose section and stag- 
nation point location are not achieved, the aft sec- 
tion of the subscale airfoil is redesigned and again 
merged with the nose section to form a new subscale 
airfoil. The flow over the new subscale airfoil is then 
analyzed and compared with that over the full scale 
airfoil The process is repeated until the desired in- 
viscid velocity distribution over the nose section and 
the stagnation point location are achieved. 

In the next step, the subscale airfoil circulation, 
water droplet trajectories and water droplet im- 
pingement characteristics are determined from AIR- 
DROP. The airfoil droplet impingement code, AIR- 
DROP, written by Bragg 12 predicts droplet trajecto- 
ries and the resultant impingement efficiency on sin- 
gle element airfoils in incompressible flow. The code 
has been validated against NACA airfoil droplet im- 
pingement data and compares well when the cloud 
droplet size distribution is modeled correctly and 
the code is run matching the airfoil lift coefficient. 12 
Comparisons with predicted and measured rime ice 
accretion show good agreement. 

The numerical procedure employed by AIRDROP 
consists of two steps. First, the flowfield around the 
airfoil is determined by Woan’s method. 18 Second, 
single water droplet trajectories are calculated from 
the trajectory equation, 12 which in nondimensional 
form contains the three additional similarity param- 
eters Ru> F r and K, apart from Re and Af. Thus, 
given Ru,F ri K, the droplet initial location, and the 
airfoil geometry, single water droplet trajectories are 


determined from the trajectory equation. 12 

The individual droplet trajectories are combined 
to calculate the local impingement efficiency /?(= 
dy 0 /dS). The impingement efficiency represents the 
dimensionless mass flux of impinging droplets at a 
point on the airfoil. Here, y 0 is the initial y displace- 
ment of an impinging droplet far ahead (x 0 = — 5c/, ) 
of the airfoil, and S is the surface length of the im- 
pact location measured from the leading edge of the 
airfoil. The AIRDROP code calculates a series of 
droplet trajectories, fits a cubic spline through the 
y 0 vs 5 data points of the impinging droplets, and 
then computes the slope of the spline at a series of 
surface positions. This slope is 0 at that surface lo- 
cation. In this paper, the y 0 vs 5 plot is referred to 
as the y*-curve and the 0 vs S plot is referred to as 
the /?-curve. And the term “impingement character- 
istics” refers to both the y 0 -curve and /?-curve. 

The impingement characteristics of both the full 
scale and subscale airfoil are then compared with 
each other. If the agreement in the impingement 
characteristics is poor, the subscale airfoil is modi- 
fied and the design process is repeated again until a 
good agreement is reached. 

As will be shown later, the amount of circulation 
plays a dominant role in determining the impinge- 
ment characteristics through its impact on the flow- 
field droplet trajectories (y 0 -curve). The expression 
for the total circulation can be derived from the re- 
lation 

L = pUf = \pU 2 cC t (1) 

which yields 



Therefore, the full-scale and subscale airfoil circular 
tion is, 

Iy# = r »« = 2 C, “ ^ 3 ’ 4 ^ 

respectively, where rj is the normalised subscale air- 
foil chord length. 

Finally, in order to obtain a physically realistic 
subscale airfoil design, consideration must also be 
given to viscous and compressibility effects to de- 
termine the true merits of the design. A discussion 
of the viscous considerations is presented in a later 
section. 

IMPLEMENTATION 

To expedite the design procedure, the Eppler 
code, PROFOIL, and AIRDROP were integrated 
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Fig. 2 The Learjet 305 (GLC 305) airfoil. 
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Fig. 3 Droplet impingement efficiency predicted 
by AIRDROP for the Learjet 305 airfoil. 

into a single computer program. Then, the con- 
straints on the subscale airfoil design were defined 
in terms of the fixed nose section geometry, velocity 
distribution over the nose section, total circulation 
and the angle of attack relative to the nose section of 
a full-scale airfoil. In order to satisfy all of the above 
constraints, numerous parametric trade studies were 
performed to help identify and isolate various key in- 
dependent design variables. These independent vari- 
ables were later identified as the pitching moment 
coefficient c mo of the airfoil from which the aft sec- 
tion of the subscale airfoil is derived, the scale rj of 
the subscale airfoil, the nose droop angle 7, and the 
upper and lower surface pressure recovery locations 
x r and z r . 

To illustrate the effects of the independent design 
variables on the subscale airfoil design, the Learjet 
305 (GLC 305) airfoil, shown in Fig. 2, was selected 
as the full-scale airfoil along with the flight and icing 
conditions listed in Table 2. At these conditions for 
the GLC 305 airfoil, AIRDROP predicts a lift coeffi- 
cient Cij t = 0.736 and the circulation Tj t = 0.368. 
Figure 3 shows the corresponding /^-curves as pre- 
dicted by AIRDROP. For VMD = 20jun, AIR- 
DROP predicts the maximum limits of impingement 
as S u = 0.0076 ( x/c = 0.0019) on the upper surface 
and Si = —0.1822 (x/c = 0.1738) on the lower sur- 
face. 


Table 2 T ypical Flight and Icing conditions. 

Airspeed, U =87 m/s (175 kt) 

Static temperature, T = -5 deg C 

Reynolds number, Re = 6 x 10 6 

Mach number, M = 0.28 

VMD = 15-40/im 

Angle of attack, a = 6 deg 

Since the limits of impingement define the surface 
within which ice will accrete on the airfoil, only that 
part of the full scale geometry need be fixed as the 
nose section for the subscale airfoil. The nose section 
size is kept to a minimum, thereby, allowing more 
flexibility in the design of the aft section to satisfy 
the constraints. Thus, the nose section geometry was 
selected as the full scale airfoil surface from x/c = 
0.05 on the upper surface to x/c = 0.20 on the lower 
surface. Moreover, a half-scale (77 = 0.5) subscale 
model was selected as the baseline case. Based on 
the size of droplets under consideration, the effect 
of gravity on the droplets was considered negligible 
and, therefore, was ignored. 

Most of the important effects can be examined by 
only considering inviscid effects; that is, boundary- 
layer displacement effects are second order relative 
to the effects of pitching moment, subscale airfoil 
chord length and the nose droop. Thus the remain- 
der of this section is divided into inviscid and viscous 
considerations. 

Inviscid Considerations 
Effect of Pitching Moment Coefficient (cm*) 

The effect of the pitching moment coefficient is il- 
lustrated in Fig. 4, in which, an increase in the pitch- 
ing moment coefficient c m0 (more negative) results 
in a subscale airfoil with a greater aft camber and, 
therefore, a higher aft loading as well as an increase 
in the amount of circulation. The change in circu- 
lation with c mo is found to be nearly linear. The 
droplet impingement characteristics, specifically the 
y*-curves, also indicate a strong dependence on the 
value of circulation which makes the pitching mo- 
ment coefficient c mo the main independent design 
variable (Note that, Cm 0 assumes a role of a depen- 
dent design variable in the design of the aft sec- 
tion). Figure 4(d) indicates that the subscale air- 
foil requires slightly less circulation (by 4.5%) than 
the full scale airfoil to achieve full-scale droplet im- 
pingement. One explanation is that the subscale air- 
foil is able to achieve full-scale droplet impingement 
with slightly less circulation due to the distribution 
of vorticity. In the case of a subscale airfoil the vor- 
ticity is more “concentrated” near the leading-edge 
than in the case of the full scale airfoil resulting in 
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Fig. 4 The effect of pitching moment coefficient 
on (a) the velocity distribution, (b) initial dis- 
placement, (c) droplet impingement efficiency 
and (d) the tangent droplet trajectories. 


a greater upwash in close proximity of the airfoil. 
Thus a lower value of overall circulation is required 
to simulate full scale droplet impingement. 

Effect of Chord Length (r/) 

To examine the effect of normalized subscale chord 
length r] on the design, subscale airfoils were de- 
signed for three different values of tj, that is, 0.5, 
0.7 and 0.9. Initially, the three subscale airfoils 
were designed such that they produced the same 
amount of circulation as the full scale airfoil. Figure 



Fig. 5 The effect of chord length on (a) the ve- 
locity distribution, (b) initial displacement, (c) 
droplet impingement efficiency and (d) the tan- 
gent droplet trajectories. 

5 shows the resulting velocity distribution impinge- 
ment characteristics and the airfoil shapes. The 
results indicate that as the scale of the subscale 
model is reduced, the aft-loading on the airfoil in- 
creases significantly in order for it to produce the 
same amount of circulation. The mismatch in the 
y c -curves, Fig. 5(b), suggests that subscale models 
require less circulation to achieve full scale impinge- 
ment characteristics. Moreover, the results also sug- 
gest that the smaller the scale, the less circulation 
required to simulate full scale droplet impingement 
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Fig. 6 Results showing (a) the velocity distri- 
bution, (b) initial displacement, (c) droplet im- 
pingement efficiency and (d) the tangent droplet 
trajectories at the matched conditions. 

characteristics. The subscale airfoils shown in Fig. 6 
were designed such that the impingement character- 
istics, specifically, the y 0 -curves were matched. The 
match in y 0 -curves was achieved by designing sub- 
scale airfoils with reduced circulation as compared 
with the ones in Fig. 5. The results also indicate 
that the amount of circulation required to simulate 
full scale droplet impingement vary from (0.955F/,) 
for 77 of 0.5 to (0.983r /,) for r] of 0.9. 

Effect of Nose Droop Angle (7) 

The effect of the nose droop angle 7, shown in 
Fig. 7, becomes obvious from Fig. 8 which illustrates 



Fig. 7 The nose droop angle. 



Fig. 8 The effect of nose droop on (a) the ve- 
locity distribution, (b) initial displacement, (c) 
droplet impingement efficiency and (d) the tan- 
gent droplet trajectories. 

the usefulness of the nose droop in reducing the high 
aft-loading on airfoils. In order to keep the angle of 
attack relative to the nose section chord constant 
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for both the full scale and subscale airfoils, the sub- 
scale airfoil with a nose droop is analyzed at an ef- 
fective angle of attack a e which takes into account 
the nose droop angle. As a result, the subscale air- 
foils with nose sections droop downwards are ana- 
lyzed at higher angles of attack than those without 
the nose droop. Figure 8 shows the results of the 
subscale airfoil design with different nose droop an- 
gles for the same value of circulation as that of the 
half-scale model without the nose droop. The results 
indicate that the nose droop results in an increase 
in the camber of the subscale airfoil and, therefore, 
the subscale airfoil circulation. Moreover, the sub- 
scale airfoils with the nose drooped downwards also 
operate at higher absolute angles of attack defined 
by a e . As a result of this increase, the impingement 
characteristics show a mismatch. By decreasing the 
amount of circulation by an appropriate amount, the 
mismatch was removed as shown in Fig. 9. The re- 
duction in the value of circulation as compared with 
that for the full scale airfoil varies from (0.95517,) 
for 7 of 0 deg to (0.89217,) for 7 of -3 deg. 

Other Effects 

The upper and lower surface pressure recovery lo- 
cations x r and *x T (see Fig. 10) control to a great 
extent the shape of the airfoil near its trailing edge. 
Although, the effect of moving the pressure recovery 
locations x r and x r results in a significant amount 
of im provement in the velocity distributions and ul- 
timately the viscous characteristics, the change in 
the droplet impingement characteristics is, however, 
small. 

The above study, based on inviscid considerations 
alone, illustrates the effect of different independent 
design variables on the subscale airfoil design. The 
results indicate that subscale airfoils require less cir- 
culation to simulate full scale airfoil droplet impinge- 
ment characteristics. The pitching moment coeffi- 
cient c m0 can be used effectively to achieve the de- 
sired amount of circulation. Since subscale airfoils 
tend toward high aft- loading in order to simulate the 
desired impingement characteristics, a nose droop 
can be used effectively to offset the high aft-loading 
to a large extent. The above study also reveals that 
subscale airfoils with a nose droop require even less 
circulation to achieve the desired impingement char- 
acteristics. Moreover, a subscale airfoil with a nose 
droop (downwards) must operate at a higher abso- 
lute angle of attack to simulate full scale impinge- 
ment characteristics over its nose section. Oper- 
ation at high absolute angles of attack makes the 
subscale airfoil highly susceptible to flow separation 
and, therefore, it becomes necessary to evaluate the 
performance by means of a viscous analysis of the 



Fig. 9 Results showing (a) the velocity distri- 
bution, (b) initial displacement, (c) droplet im- 
pingement efficiency and (d) the tangent droplet 
trajectories at the matched conditions. 

flowfield over the subscale airfoil at the design con- 
ditions. 

Viscous Considerations 
To determine the true merits of the design, a vis- 
cous flowfield analysis must form an essential part 
of the design. For the purpose of viscous flowfield 
analysis, XFOIL was utilized. XFOIL is a modified 
version of the ISES code 19 which has been success- 
fully applied to the design and analysis of airfoils for 
various applications varying from human-powered 
aircraft 20 to high Reynolds number transonic trans- 
port. XFOIL utilizes a fully compatible laminar and 
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Fig. 10 The upper and lower surface match and 
pressure recovery locations. 

turbulent viscous formulation, a reliable transition 
formulation and a global Newton iteration method 
to converge onto a flowfield solution. 

The viscous analysis was performed to determine 
the effect of the presence of the viscous boundary on 
the flowfield. Typically inviscid flowfield codes over 
predict the airfoil lift-curve slope and the lift on an 
airfoil as compared to viscous flowfield codes since, 
in a viscous flowfield, the presence of boundary-layer 
decambers the airfoil and, therefore, reduces the C\. 
This results in an error in the droplet trajectory cal- 
culation since at the design angle of attack, the in- 
viscid flowfield is for a higher lift coefficient, and 
therefore, greater circulation. To account for this 
effect, a procedure called the “Matched Lift Coef- 
ficient Method” is employed, in which, the inviscid 
flowfield is analyzed at matched lift coefficient in- 
stead of matched angle of attack with the viscous 
flowfield. A brief outline of this procedure as ap- 
plied to the subscale airfoil design is as follows. 

Initially, the viscous Cij t is determined at the de- 
sign angle of attack with the help of XFOIL. Using 
Cjj,, an angle of attack <*ij$ is found such that 
running the inviscid flowfield code at a *,/, produces 
an inviscid (7/ which matches C/j,, the viscous C\. 
Next, the inviscid flowfield as well as the droplet 
impingement characteristics of the full scale airfoil 
are determined at oc{j s and set as the target for 
the subscale airfoil design. A subscale airfoil is then 
designed to match the target flowfield and impinge- 
ment characteristics. Once a match is achieved, a 
viscous analysis of the subscale airfoil is performed 
at the matched conditions to determine the viscous 
Cf,„. As in the full scale airfoil case, an inviscid 
a,- , , is calculated and is used to determine the invis- 
cid flowfield and droplet impingement characteristics 
for comparison with the target flowfield and droplet 
impingement characteristics. If the desired charac- 


Table 3 Design Flight and Icing conditions. 

Variable 

Full scale 

Subscale 

U, m/s 

87 

87 

r, deg c 

-5 

-5 

Re 

6 x 10 6 

3 x 10® 

M 

0.28 

0.28 

c, m 

1.0 

0.5 

VMD, fim 

20 

20 

a, deg 

6 

6 

T. de 8 

0 

-3 

a e , deg 

6 

9 

Table 4 

The converged solution. 


Dependent Variables Independent Variables 


Ks = 

0.3 

<t>u 

= 189.53 deg 

Crno = 

: -0.065 

Vl 

= 2.133 

y(x m 

= 0.05) 

a* 

= 8.93,11.93,14.93 deg 

y(x m 

= 0.20) 

cT 

= 1.17 deg (all segments) 

r = 6 

deg 


= 0.0114c x r = 0.4746c 


teristics are achieved, the design is complete, other- 
wise, the subscale airfoil is modified and the whole 
process is repeated again until the desired match is 
achieved. 

A DESIGN EXAMPLE 

In this section, a specific design example is pre- 
sented with the objective to design a half-scale model 
of the GLC 305 airfoil that simulates full scale 
droplet impingement. Table 3 lists the flight and 
icing conditions for the final design, whereas, Ta- 
ble 4 lists the final values of the design variables for 
the converged solution. The subscale airfoil was de- 
signed with a finite trailing-edge angle r = 6 deg. 
The effects due to compressibility were also consid- 
ered during the viscous flow analysis of both the air- 
foils. 

Fig. 11 shows the comparison between the in- 
viscid velocity distributions for the converged solu- 
tion at = 6 deg. Fig. 12 shows the compari- 
son between the velocity distributions (viscous) at 
the design conditions, where, C\j 9 = 0.7690 and 
C l$9 - 0.6074. The respective inviscid velocity dis- 
tributions for the matched lift coefficient case are 
shown in Fig. 13(a). All the figures show good agree- 
ment in velocity distribution over the common nose 
section. The comparison of the impingement char- 
acteristics corresponding to the respective matched 
lift coefficient cases is shown in Fig. 13(b) and (c), 
whereas, a comparison of tangent droplet trajecto- 
ries is shown in Fig. 13(d). The results indicate ex- 
cellent agreement in impingement efficiency. The 
tangent droplet trajectories, although originating at 
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x/c 


Fig. 11 Comparison between the inviscid velocity 
distribution at ct e = 6 deg. 



-0.2 0 0.2 0.4 0.6 0.8 1 1.2 


x/c 

Fig. 12 Comparison between the velocity distri- 
bution at the design conditions listed in Table 3. 

different locations along the y-axis, are matched in 
the vicinity of the leading-edge. This is consistent 
with the observations made during the case studies 
that subscale airfoils require a lower value of cir- 
culation to achieve full scale droplet impingement 
characteristics. 

CONCLUSIONS 

Several important conclusions can be drawn from 
this study. First, it is shown that subscale air- 
foils with full-scale leading edges can be designed 
to exhibit full-scale droplet impingement and, there- 



Fig. 13 Comparison between (a) the velocity 
distribution, (b) the initial displacement, (c) 
droplet impingement efficiency and (d) the tan- 
gent droplet trajectories at <*«,/# and cor- 

responding to the respective matched lift coeffi- 
cients. 

fore, ice accretion. Second, the results indicate that 
subscale airfoils require less circulation to simulate 
full scale airfoil droplet impingement characteristics. 
The pitching moment coefficient, of the airfoil from 
which the aft section for the subscale airfoil is de- 
rived, can be used effectively to achieve the desired 
amount of circulation on the subscale airfoil. Third, 
since subscale airfoils tend toward high aft-loading 
in order to simulate the desired droplet impingement 
characteristics, a nose droop can be used effectively 
to offset the high aft-loading. Fourth, an airfoil with 
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Full Scale 
Subscale 



Fig. 14 The Anal subscale airfoil and the Learjet 
305 airfoil. 

a nose droop (downwards) must operate at a higher 
absolute angle of attack in order to keep the same 
angle of attack relative to its nose section as the full 
scale airfoil to simulate full scale impingement char- 
acteristics. Operation at high absolute angles of at- 
tack makes the subscale airfoil highly susceptible to 
flow separation and, therefore, it becomes necessary 
to integrate the viscous analysis of the flow field over 
the subscale airfoil into the design process. Fifth, to 
incorporate viscous and compressibility effects, the 
“matched lift coefficient method” outlined in the pa- 
per was applied successfully in the final design ex- 
ample. 

Although, the design method outlined in this pa- 
per is only limited to a point design, the method can 
be extented to a multipoint design similar in lines to 
the existing multipoint inverse airfoil design meth- 
ods by integrating viscous boundary-layer equations 
and the droplet trajectory equation with the inverse 
airfoil design method. 
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ABSTRACT 

A design procedure for hybrid airfoils with full- 
scale leading edges and redesigned aft-sections that 
exhibit full-scale airfoil water droplet impingement 
characteristics throughout a given C \- range is pre- 
sented. The design procedure is an extension of the 
method first published by Saeed, et al., in that it 
not only allows for subcritical and viscous flow anal- 
ysis in the design but is also capable of off-design 
droplet impingement simulation through the use of 
a flap system. The limitations of the flap-system 
based design for simulating both on- and off-design 
full-scale droplet impingement characteristics and 
surface velocity distribution are discussed with the 
help of specific design examples. In particular, the 
paper presents the design of two hybrid airfoils at 
two different angles of attack, such that they simu- 
late both full-scale velocity distribution as well as 
droplet impingement at the respective design an- 
gles of attack. Both of the hybrid airfoils are half- 
scale airfoil models with a 5% upper and 20% lower 
full-scale surface of the Learjet 305 airfoil leading- 
edge. The effect of flap deflection and droplet size 
on droplet impingement characteristics is also pre- 
sented to highlight the important limitations of the 
present method both on and off design. The paper 
also discusses important compromises that must be 
made in order to achieve full-scale ice accretion sim- 
ulation throughout a desired C \ -range and suggests 
alternatives such as applying a multipoint design ap- 
proach for the design. 

NOMENCLATURE 

Ci = airfoil lift coefficient 

c = airfoil chord length 

5 = airfoil surface arc length measured from 

the leading-edge 

T = freestream static temperature 
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V = surface velocity 

Koo = freestream velocity 

V = surface velocity normalized by 

VMD = volume median droplet diameter 
x, y = airfoil coordinates 

a = angle of attack relative to the chord line 
a c = effective angle of attack relative to the 
nose-section chord line, a — 7 

0 = local impingement efficiency 

7 = nose droop angle 

T = circulation strength normalized by V^c 

6 = droplet diameter 

5f = flap deflection, deg 

Subscripts; 

fs — full-scale airfoil 

1 = lower surface 

ss = subscale airfoil 

u = upper surface 

INTRODUCTION 

Recent aircraft accidents have raised important 
flight safety 1 ' 5 issues related to the operation of air- 
craft under severe weather conditions. To improve 
flight safety, a better understanding of the effect 
of ice accretion on the aerodynamic performance of 
modern airfoils is required. One important step in 
the process is to evaluate the aerodynamic perfor- 
mance of the airfoil sections, or the wing as a whole, 
at the icing conditions within the certification ic- 
ing envelop that results in the largest performance 
penalties. 

For aircraft safety, one of the most important 
performance parameters is the maximum lift coef- 
ficient. Therefore, while drag and pitching moment 
are important, the icing condition that results in the 
largest degradation in maximum lift coefficient is the 
most critical icing condition. The determination of 
the critical ice accretion and its aerodynamic effect 
on a set of modem airfoils, typical of those in use 
on aircraft, is underway at NASA Lewis Research 
Center. The research reported here is part of this 
larger effort. 

Owing to the difficulties and uncertainties in ice 
accretion scaling , 6 ' 14 testing at full-scale is desir- 
able, yet costly. Moreover, available ice accretion 
tunnels are too small to test full-scale airfoils or 
wings of most aircraft of interest. One way to ex- 


1 



pand the usefulness of existing icing tunnels and to 
facilitate testing of aircraft deicing/anti-icing sys- 
tems is to test “hybrid airfoils” or “sub-scale air- 
foils” with full-scale leading edges and redesigned 
aft sections to provide full-scale icing conditions at 
the leading edge. The term “hybrid method” refers 
to using a full-scale leading edge to match the full- 
scale ice accretion. The aft section of the hybrid 
airfoil is specially designed to provide flowfield and 
droplet impingement similar to that on the full-scale 
airfoil leading-edge. One such approach 15 used air- 
foils with full-scale leading edges and truncated aft- 
sections to simulate the flowfield of the full scale, 
thereby avoiding the ice-accretion process on the air- 
foil leading edge and the associated seeding issues 
altogether. Interestingly, neither the approach, nor 
its range of application, received much attention de- 
spite its numerous merits since it permits an indepth 
study of droplet impingement and ice accretion on 
full-scale leading-edge sections within the capabili- 
ties of current icing research facilities. 

In the absence of a systematic study to provide 
insight into the design of the aft section, a recent 
study 16 was carried out in which a design proce- 
dure for hybrid airfoils was successfully developed 
and demonstrated with specific design examples. 
The study showed that hybrid airfoils could be de- 
signed to exhibit both the full-scale velocity distribu- 
tion on its nose section as well as full-scale droplet- 
impingement characteristics and, therefore, ice ac- 
cretion. An inherent limitation of the design pro- 
cedure outlined in the study 16 is that the method 
was restricted to a single-point design and, there- 
fore, lacked the capability to handle off-design cases. 
Moreover, the method used the “matched lift coeffi- 
cient” technique to correct for viscous effects. 

To overcome these limitations, the present study 
was carried out with the objective to expand the 
scope of the single- point design procedure of Ref. 16 
to a method that enables the hybrid airfoils to ex- 
hibit both full-scale velocity distribution as well as 
droplet impingement and, therefore, ice accretion 
throughout a desired C \— range or a range of angles 
of attack a. 

The task of simulating off-design full-scale droplet 
impingement, as will be shown later, is successfully 
accomplished by introducing a plain flap on the hy- 
brid airfoil. The use of a plain flap, however, fails 
to simulate full-scale velocity distribution at the off- 
design conditions. Since the difference in the veloc- 
ity distribution on the nose section will effect the 
thermodynamics of ice accretion as the droplets im- 
pinge on the surface, it therefore becomes necessary 
to simulate the full-scale velocity distribution in ad- 


dition to droplet impingement at the off-design con- 
ditions. Thus, to simulate both the full-scale ve- 
locity distribution as well as droplet impingement 
on the nose section of the hybrid airfoil throughout 
a desired C\ —range, it is necessary to formulate a 
multipoint hybrid airfoil design method. 

To set the stage for the multipoint design method, 
the paper presents the design of two half-scale hy- 
brid airfoils that are designed at two different angles 
of attack such that they simulate both the full-scale 
velocity distribution as well as droplet impingement 
characteristics on the nose sections at their respec- 
tive design angles of attack. The velocity distribu- 
tion and droplet impingement characteristics of the 
two hybrid airfoils are then analyzed at an off-design 
angle of attack and compared with that of the full- 
scale airfoil. The results are then used to highlight 
the limitations of the present method and, therefore, 
suggest a need for a multipoint design method. Im- 
portant compromises that must be made to achieve 
a multipoint design for full-scale ice accretion sim- 
ulation throughout a desired C \ -range are also dis- 
cussed. 

DESIGN APPROACH 

The hybrid airfoil design procedure for full- 
scale flowfield and droplet impingement simulation 
uses validated computational airfoil aerodynamics 
and droplet impingement codes, 17-32 specifically, 
an inverse design method, 32 the Eppler code, 25-27 
XFOIL 31 and AIRDROP. 20-22 Reference 16 gives a 
brief discussion on each of these codes. For a more 
detailed discussion, the reader is referred to the as- 
sociated literature. 

Unlike the method presented in Ref. 16 where 
in the potential flow is corrected for viscous ef- 
fects using the “matched lift coefficient” technique, 
the present method uses a modified version of 
XFOIL. The modified version of XFOIL was ob- 
tained by integrating the droplet-trajectory and 
impingement-characteristics calculation subroutines 
from the AIRDROP code into the XFOIL code. This 
was especially done to take advantage of XFOIL’s 
ability to analyze both inviscid/viscous flow as well 
as incompressible/subcritical flows unlike the AIR- 
DROP code, which is purely based on incompress- 
ible flow formulation. In this paper, the modi- 
fied version of XFOIL will be referred to as the 
XFOIL/ AIRDROP code. Once the flowfield is de- 
termined using known flight and icing conditions, 
the droplet trajectory calculation subroutines are 
then used in conjunction with the flow solver subrou- 
tines to determine the water droplet impingement on 
the airfoil surface. 
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A conceptual illustration of the hybrid airfoil de- 
sign procedure is shown in Fig. 1. A brief summary 
of these steps is as follows. First, a full-scale air- 
foil geometry is selected and the desired flight and 
icing conditions are specified. In particular, the 
Learjet 305 airfoil (shown in Fig. 2) is used in this 
study to demonstrate the design procedure. The 
XFOIL/ AIRDROP code is then used to predict the 
limits of droplet impingement. Once the limits of im- 
pingement are known on the leading edge of the full- 
scale airfoil, that part of the full-scale airfoil geome- 
try is fixed for the subsequent hybrid airfoil shapes. 
As in Ref. 16, this fixed leading-edge section will be 
referred to as the nose section and the remaining 
section of the subscale airfoil profile will be referred 
to as the aft section. The aft section of the hybrid 
airfoil is then designed to provide full-scale flowfield 
and droplet impingement characteristics on the nose 
section of the hybrid airfoil. 

An initial geometry for the aft section is obtained 
through the use of a multipoint inverse airfoil design 
code 32 (PROFOIL). The design of the intermediate 
airfoil, from which the aft section of the subscale 
airfoil is derived, is governed by several constraints, 
namely, the scale of the subscale airfoil, the upper 
and lower surface thickness and slope at the junction 
between the nose and aft sections, and a desired form 
for the pressure recovery characteristics. Apart from 
these constraints, additional continuity and closure 
constraints that form an integral part of the inverse 
design methodology 32 are also satisfied to achieve 
a physically realizable design. A multi-dimensional 
Newton iteration scheme is further employed to sat- 
isfy these constraints. 

The flow over the hybrid airfoils is then analyzed 
using the XFOIL/AIRDROP code. In order to have 
a physically similar flow in the vicinity of the nose 
section of both the hybrid and the full-scale airfoils, 
the analysis is performed at the same angle of attack 
relative to the nose-section chord of both the air- 
foils. The local velocity distributions over the nose 
section and the stagnation point locations on both 
the hybrid and full-scale airfoils are then compared. 
If the desired velocity distribution over the nose sec- 
tion and stagnation point location are not achieved, 
the aft section of the hybrid airfoil is redesigned and 
again merged with the nose section to form a new 
hybrid airfoil. The flow over the new hybrid airfoil 
is then analyzed and compared with that over the 
full-scale airfoil. The process is repeated until the 
desired velocity distribution over the nose section is 
achieved. 

In the next step, the water droplet trajectories 
and water droplet impingement characteristics are 


determined from the XFOIL/AIRDROP code. The 
individual droplet trajectories are combined to cal- 
culate the droplet impingement characteristics of the 
airfoil. The droplet impingement characteristics of 
both the full-scale and the hybrid airfoil are then 
compared. If the agreement in the droplet impinge- 
ment characteristics is poor, the hybrid airfoil is 
modified and the design process is repeated again 
until good agreement is reached. At this stage, the 
single-point design is accomplished. To achieve off- 
design full-scale ice accretions or droplet impinge- 
ment characteristics, a plain flap is employed on 
the hybrid airfoil. Thus, by deflecting the flap, 
the desired droplet impingement characteristics are 
achieved over the hybrid airfoil for the off-design 
cases. 

The off-design cases reveal, as will be shown in the 
next section, certain important limitations of the de- 
sign method. These limitations include 1) the onset 
of flow separation on the hybrid airfoils at moder- 
ate to high angles of attack conditions and 2) a mis- 
match in the velocity distribution on the nose section 
at off-design angles of attack. The former limitation 
can be improved either by using a more sophisti- 
cated flap system or by applying less conventional 
techniques such as boundary-layer control through 
slot suction 33,34 or circulation control via trailing- 
edge blowing. The latter, however, is an important 
limitation of the present design method and can be 
overcome by using a multipoint design approach. 

IMPLEMENTATION 

In this section, the effects of various parameters on 
two single-point airfoil designs are discussed. In par- 
ticular, two half-scale hybrid airfoils were designed 
at different angles of attack such that they simu- 
lated both the full-scale velocity distribution on the 
nose section as well as droplet impingement char- 
acteristics at the design conditions (single-point de- 
sign). The off-design full-scale velocity distribution 
and droplet impingement simulation characteristics 
of each hybrid airfoil are compared to highlight im- 
portant limitations of the present method. 

Single-Point Design and Simulation 

The design of two half-scale models of the GLC 
305 airfoil that simulate full-scale velocity distribu- 
tion and droplet impingement is presented. Of the 
two hybrid airfoils A and B, hybrid airfoil A is de- 
signed to simulate full-scale ice accretion at a = 2 
deg while hybrid airfoil B is designed to simulate 
full-scale ice accretion at a = 6 deg along with the 
icing conditions: Voo = 90 m/s (175 kt), T = — 10°C, 
Re = 6 x 10 6 , M = 0.28 and VMD = 20/xm. While 
it is realized that in flight the conditions will change 
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Table 1 Design flight and icing conditions. 


Variable 

Full scale 

Hybrid A 

Hybrid B 

Voo, rn/s 

90 

90 

90 

T, deg C 

-10 

-10 

-10 

Re 

6 x 10 6 

3 x 10 6 

3 x 10 6 

M 

0.28 

0.28 

0.28 

c, m 

1.0 

0.5 

0.5 

VMD , /xm 

20 

20 

20 

a, deg 

2,6 

2 

6 

7. deg 

0 

-1.5 

-3 

a«, deg 

2,6 

3.5 
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with angle of attack, the conditions for both an- 
gles are held constant here to simply illustrate the 
method. 

As a first step, the droplet impingement efficiency 
0 for the GLC 305 airfoil corresponding to the 
given flight and icing conditions is determined by 
the XFOIL/AIRDROP code. The results are shown 
in Fig. 3. For a = 6 deg, the XFOIL/AIRDROP 
code predicts the maximum limits of impingement 
as S u — 0.0076 (x/c = 0.0019) on the upper sur- 
face and Si = -0.1822 (x/c = 0.1738) on the lower 
surface. Since the limits of impingement define the 
surface over which ice will accrete on the airfoil, only 
that part of the full-scale airfoil geometry needs be 
fixed as the nose section for the hybrid airfoil. Thus, 
the nose-section geometry for both the hybrid air- 
foils was selected as the full-scale airfoil surface from 
x/c = 0.05 on the upper surface to x/c = 0.20 on the 
lower surface. The two hybrid airfoils were then de- 
signed following the procedure illustrated in Fig. 1. 
Table 1 lists the flight and icing conditions for the 
final single- point design. 

A comparison of the full-scale airfoil velocity dis- 
tribution with that of the individual hybrid airfoil 
velocity distributions (Figs. 4a and 5a) at the single- 
point design conditions shows good agreement over 
the common nose section. Comparisons of the im- 
pingement characteristics (Figs. 4b and 5b) and tan- 
gent droplet trajectories (Figs. 4c and 5c) also indi- 
cate excellent agreement with that of the full-scale. 
The tangent droplet trajectories, although originat- 
ing from different locations upstream are matched in 
the vicinity of the leading edge. This is consistent 
with the observations made during the case studies 
in Ref. 16. At this point, the single-point design for 
full-scale velocity distribution and droplet impinge- 
ment simulation is complete and the two hybrid air- 
foils along with the Learjet 305 airfoil are shown in 
Fig. 6. 

Effect of Droplet Size 

The impingement characteristics, i.e., the lim- 


its of impingement, the impingement efficiency 0 
(0— curve) and the maximum point on the 0— curve, 
referred to as 0max > of an airfoil depend to a large 
extent on the size of the water droplets in the flow. 
In the case of small droplets, the droplet drag domi- 
nates and the particle is very responsive to the flow- 
field acting almost as a flow tracer; whereas, in the 
case of large droplets, the droplet inertia dominates 
and the particle is less sensitive to changes in the 
flowfield. Thus, an increase in the droplet size re- 
sults in an increase in the impingement efficiency /?, 
0max and the limits of impingement. It, therefore, 
becomes necessary to examine the effect of different 
droplet size on full-scale droplet impingement sim- 
ulation. Since, in an actual icing cloud, the water 
droplets have diameters ranging from 5-50 /im, the 
impingement characteristics of the hybrid airfoil A 
were determined for two different droplet sizes. The 
results are presented in Fig. 7 and show good agree- 
ment where the droplet sizes are less than that se- 
lected for the single point design. For larger droplet 
size, a good overall agreement can be seen, however, 
the limits of impingement and 0 ma x differ slightly. 

Off-Design Simulation 

To simulate full-scale ice accretion or droplet im- 
pingement characteristics throughout a desired Ci — 
or a— range, a flap system was employed on each of 
the hybrid airfoils. The objective was to match both 
the velocity distribution as well as the droplet im- 
pingement characteristics at any off-design angle of 
attack by an appropriate amount of flap deflection. 
To accomplish this task, the two hybrid airfoils were 
analyzed at off-design angles of attack, in partcular, 
the hybrid airfoil A designed to simulate conditions 
at a = 2 deg was analyzed at a = 6 deg while the 
hybrid airfoil B designed to simulate conditions at 
a = 6 deg was analyzed at a = 2 deg. The results 
are shown in Figs. 8 and 9 in which the hybrid airfoil 
velocity distribution and impingement characteris- 
tics are shown with and without the appropriate flap 
deflection necessary to simulate full-scale droplet im- 
pingement. The results show that, although the use 
of a flap on hybrid airfoils can be very effective in 
simulating full-scale droplet impingement character- 
istics at an off-design condition, it is, however, not 
able to accurately simulate full-scale velocity distri- 
bution over the nose section of that hybrid airfoil. 

To determine the optimum flap setting, the root- 
mean-squares difference in local impingement effi- 
ciency RMSp and in normalized surface velocity 
RMS V were calculated for different angle of at- 
tack and flap settings. Mathematically, RMSp and 
RMS V are defined as 
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RMS 0 = \\p f ,(S)-0„(S)\\ (1) 

RMS v =\\V Ia (S)-V >t (S)\\ (2) 

where Si < S < S u . 

Figures 10a, b show the variation in RMSp and 
RMSy ■, respectively, for different angles of attack 
and flap settings 6 / for the hybrid airfoil A designed 
for a = 2 deg while Figs. lla,b show similar plots 
for the hybrid airfoil B designed for a = 6 deg. The 
optimum flap deflection was then selected as the one 
that corresponds to the minimum value of RMSp. 

The optimum flap settings corresponding to each 
angle-of-attack case are plotted in Fig. 12a for clar- 
ity. Figure 12b, on the other hand, shows a com- 
parison of the circulation T of both the hybrid air- 
foils with that of the full-scale. The results indicate 
that the hybrid airfoils require less circulation than 
the full-scale airfoil to simulate full-scale droplet im- 
pingement and that the difference between the full- 
scale and hybrid airfoil circulation is nearly constant 
until significant flow separation occurs on the hybrid 
airfoils. Beyond this point, the hybrid airfoil circu- 
lation starts to fall off gradually and, therefore, sug- 
gests the limit to which a hybrid airfoil can be used 
to simulate full-scale droplet impingement charac- 
teristics. 

It is important to note in Figs. 10 and 11 that the 
RMS^r values are an order of magnitude higher than 
the corresponding RMS@. Although contributions 
to the RMS values due to numerical noise cannot 
be ruled out completely, differences in surface veloc- 
ity will certainly effect the thermodynamics of ice 
accretion. Thus, it becomes necessary to incorpo- 
rate the ice accretion process in the design method 
in addition to flow and droplet impingement analy- 
sis. 

The effect of larger droplet size on off-design sim- 
ulation is shown in Fig. 13. Similar trends can be 
observed as in the on-design case. Since large sized 
droplets result in an increase in the limits of im- 
pingement, they together with the angle of attack of 
interest may dictate the size of the nose section and, 
thus, limit the range of application of the present 
method. 

CONCLUSIONS 

Several important conclusions can be drawn from 
this study. First, it is shown that it is possible to de- 
sign hybrid airfoils with full-scale leading edges and 
redesigned aft-sections that exhibit full-scale airfoil 
water droplet impingement characteristics through- 
out a given Q— range. The results indicate the use- 
fulness of a flap system in simulating off-design full- 
scale droplet impingement characteristics. The use 


of flap for full-scale droplet impingement simulation 
is, however, restricted to low and moderate angles 
of attack since at high absolute angles of attack to- 
gether with high flap deflections, the hybrid airfoils 
become susceptible to flow separation. This limita- 
tion can, however, be overcome by the use of a more 
sophisticated flap system or by the application of 
boundary-layer control methods. 

The results of off-design simulation also reveal the 
existance of small differences in surface velocity dis- 
tribution within the limits of droplet impingement. 
Since this difference in surface velocity will affect the 
thermodynamics of ice accretion and prevent full- 
scale ice accretion simulation, the present method 
should be modified to include the effects of ice ac- 
cretion as well into the design of hybrid airfoils. 
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(a) 



Fig. 1 The subscale/hybrid airfoil design 
procedure. 



Fig. 2 Lear jet 305 (GLC 305) airfoil. 
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Fig. 3 Droplet impingement efficiency for the 
Learjet 305 airfoil. 



Fig. 4 Comparison of full-scale and hybrid airfoil A 
(a) velocity distributions, (b) droplet impingement 
efficiencies and (c) tangent droplet trajectories at 
the design angle of attack a = 2 deg. 
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Fig. 5 Comparison of full-scale and hybrid airfoil B 
(a) velocity distributions, (b) droplet impingement 
efficiencies and (c) tangent droplet trajectories at 
the design angle of attack a = 6 deg. 


Full-scale 

Hybrid A 

Hybrid B 



Fig. 6 The two hybrid airfoils and the Learjet 305 
airfoil. 



Fig. 7 Comparison of full-scale and hybrid airfoil A 
droplet impingement efficiencies for (a) 5 micron and 
(b) 40 micron droplet size at the design angle of attack. 
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Fig. 8 Comparison of full-scale and hybrid airfoil A 
(a) velocity distributions and (b) droplet impingement 
efliciencies at ofT-design angle of attack a = 6 deg with 
and without flap deflection. 


Fig. 9 Comparison of full-scale and hybrid airfoil B 
(a) velocity distributions and (b) droplet impingement 
efficiencies at off-design angle of attack a=* 2 deg with 
and without flap deflection. 
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Fig. 10 The variation in the RMS values for 
different angles of attack and flap settings 
for the hybrid airfoil A. 



Fig. 11 The variation in the RMS values for 
different angles of attack and flap settings 
for the hybrid airfoil B. 
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Fig. 12 Plot of (a) the optimum flap deflection 
and (b) the respective airfoil circulation. 
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Fig. 13 Comparison of off-design droplet impinge- 
ment efficiency for 40 micron droplet size at ct = 4 deg. 
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